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ABSTRACT 

We study the H 2 molecular content in high redshift damped Lyman-a systems (DLAs) as a function of the H i column density. We find 
a significant increase of the H 2 molecular content around log A(Hi) (cm“^) ~ 21.5 - 22, a regime unprobed until now in intervening 
DLAs, beyond which the majority of systems have log A(H 2 ) > 17. This is in contrast with lines of sight towards nearby stars, where 
such H 2 column densities are always detected as soon as log N(H i) > 20.7. This can qualitatively be explained by the lower average 
metallicity and possibly higher surrounding UV radiation in DLAs. However, unlike in the Milky Way, the overall molecular fractions 
remain modest, showing that even at a large N(H i) only a small fraction of overall H i is actually associated with the self-shielded 
H 2 gas. Damped Lyman-a systems with very high-A(H i) probably arise along quasar lines of sight passing closer to the centre of 
the host galaxy where the gas pressure is higher. We show that the colour changes induced on the background quasar by continuum 
(dust) and line absorption (Hi Lyman and H 2 Lyman & Werner bands) in DLAs with log A(Hi) ~ 22 and metallicity ~1/10 solar is 
significant, but not responsible for the long-discussed lack of such systems in optically selected samples. Instead, these systems are 
likely to be found towards intrinsically fainter quasars that dominate the quasar luminosity function. Colour biasing should in turn be 
severe at higher metallicities. 
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1. Introduction 


The atomic to molecular hydrogen transition is a prerequisite 
process for star formation through the collapse of molecular 
clouds and there fore has important implica tions for the evolution 
of galaxies (e.g. lKennicutt & Evansll201^ . The relative amount 
of dense molecular and diffuse atomic gas in nearby galaxies is 
found to be cor related with the hydrostatic pressure at the galac¬ 
tic mid-plane (iBlitz & Rosolowskvl l2006l) . which is driven by 
the gravity of gas and stars. This is a natural consequence of 
thermal equilibrium of the gas, leading to mult iple phases un¬ 
der an external pressure (e.g. lWolfire et al.l|l995l) . The transition 
between H i and H2 can then be linked to a critical gas surface 
mass density above which star f ormation is triggered, inducing a 
Schmidt-Kennicutt relation (e.g. Schavel2001tlAltav et al.ll201 ll: 


Lagos et al.ll201 lUPonning et al.ll2014t)y 


The local abundance of H2 in the interstellar medium (ISM) 
depends on the balance betwe en its form ation, primarily on 
the surface of dust grains (e.g. I.Iuralll974 but also in the ga s 
phase though the H -hH—»H 2-i-e“ reaction. iBlack et'^ll987h . 
and its dissociation by UV photons. Because the dissociation 
occurs through Lyman and Werner band line transitions (e.g. 
[Palgarno & Stephenslll97(ih . self-shielding becomes very effi¬ 
cient when H2 absorpti on lines from several ro tational levels 
become saturated (e.g. iDraine & Bertoldil Il996h . Dust grains 
also absorb Lyman and Werner band photons further contribut¬ 
ing to decreasing the photo-dissociation rate. Theoretical mi¬ 
crophysics models that include detailed treatment of the forma¬ 
tion of H2 onto dust grains and the dust- and self-shielding of 
H2 show that the conversion from atomic to molecular occurs 
above a V( H i)-threshold that incr e ases with decreasing metal- 
licity (e.g. Krumholz et al. 2009I: M cKee & &umhola I2010I: 


iGnedin & &avtsovll2Ql 1: Sternberg et al.ll20~ ). 


A sharp increase in the H2 column densities has been first 
noticed above log V (H i) ~ 20.7 in the local Galactic ISM by 
ISavage et al.l (Il977[). In turn, the first studies of the Magellanic 
clouds bv lTumlinson et al.l (l2002h did not reveal any dependence 
of the H2 content on the H i column density. This was explained 
by a high average UV radiation due to intense local star forma¬ 
tion activity together with lower metallicities. At high redshift, 
H2 is generally detect ed in about 10% of damped Lyman-o' sys¬ 
tems (DLAs) or less (|Petitiean et al.ll2000l: iLedoux et al.ll2003l 
iNoterdaeme et al.ll2008h . Physical conditions in these sub-solar 
metallicity systems indicate densities of the order of n ~ 50 cm~^ 
in the cold neutral mediu m, and ambient radiation field a few 
times t he Draine field (e.g .JSrian and et al.ll2005HNeeleman et~an 
I2OI5I) . INoterdaeme et al.l (l2008l) noted that the presence of H2 
does not strongly depend on the total neutral hydrogen col¬ 
umn density up to logV(Hi) ~ 21.5. They concluded that large 
molecular hydrogen content, as predicted bv lSchav3(l2001[) may 
be found at higher column densities. 

Here, we investigate the H2 content of high redshift Hi- 
selected DLAs at the extreme H i column density end, a regime 
almost unprobed until now and only made reachable very re¬ 
cently thanks to very large DLA datasets (INoterdaeme et al.l 
I2OI2I) and high-resolution spectroscopic follow-up. 


2. The atomic to molecular hydrogen transition 

We have recently searched for H2 in fou r extremely strong DLAs 
(ESDLAs, defined as log V(Hi) > 21.7. INoterdaeme et al.l2014l) 
using the Ultraviolet and Visual Echelle Spectrograph (UVES) 
on the Very Large Telescope (VLT). This brings the number of 
ESDLAs with H2 searches (all with VLT/UVES) to seven. De¬ 
tails of Hi and H2 measurements in these ESDLAs are sum¬ 
marised in Table[T] Combining this with other measurements we 
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Fig-1- Left: Column den sity of H 2 as a function of that of H i. Data for high- redshift (z > 1.8) quasar-DLAs (black points and arrows) are from 


Noterdaeme et al.l j2008L updated to more recent determi nation in a few case s . Ilvanchik et alj2 01(]|:lRahmani et aLl l2013l ; |Albornoz Vasguez et al.l 
20141) . A dditional measureme nts in ESDL As are from iNoterdaerne et ahl ll201.5h arid' Guim^ese^^.1 ll20121) 7 Data for G RB-DLAs ( green) 
are from iLedoux et all ) l2009l upper limits). IProchaska et alj ) l2009l GRB080607). iKriihler etaLn2oTA GRB120815A), and iFriis et alj fcoiSL 
GRB121024A). With the exception of GRB080607 and GRB121024A, in which damped H 2 lines allow for accurate column density measurement 
at medium spectral resolution, all high-z values (quasar and GRB DLAs) are from UVES data. Purple (resp. orange) points represent measurements 
in the Milky Way (resp. SMC). The dotted, dashed, and dash-dotted lines represent average molecular fractions of, respectively, / = 0.1, 1, and 
10%. Right: Frequency of H 2 detection as a function of the Hi column density. For each Ai(Hi) bin, systems with log Ai(H 2 ) >17 are coded in 
colour, and systems with log Ai(H 2 ) < 17 are shown in white. Systems with upper limits that are not stringent (i.e. above log A(H 2 ) = 17) are not 
taken into account. The numbers in each box indicate the total number of systems contributing to the bin. Crossed boxes have no statistics. 


Table 1. H 2 in ESDLAs with log A(Hi)>21.7 


Quasar 

^abs 

log A(Hi) 

log (V(H 2 ) 

Ref. 

HE 0027-1836 

2.402 

21.75+0.10 

17.43+0.02 

1 

QJ0154H-1935 

2.251 

21.75+0.15 

~ 18 

2 

Q0458-0203 

2.040 

21.70+0.10 

< 14.60 

3 

QJ0816-H1446 

3.287 

22 .00+0.10 

18.66+0.30 

4 

Q 1157-1-0128 

1.944 

21.80+0.10 

< 14.50 

1 

QJ1456-hl609 

3.352 

21.70+0.10 

17.10+0.09 

2 

QJ2140-0321 

2.340 

22.40+0.10 

20.13+0.07 

2 


Re ferences, ('ll [ R ahtnani et alJ ([20]j[hj(2) jNoterdaeme et ^ (1201 .51) : 
(3) lNoterdaeme eTaiTil^OOsI) : (4) [Guimaraes 61^11)1^121) . 


explore the H 2 content as a function of Ai(H i) in DLAs while 
refraining from drawing any conclusion on the overall H 2 detec¬ 
tion rate. We can do so since the DLAs used for this study were 
selected only on the basis of their neutral hydrogen content. For 
this reason, we do not include recent H 2 detections obtained by 
directly targeting systems based on the presence of cold ga£]. 

In the left panel of Fig. [1] we compare the total H 2 column 
dens ity versus that of Hi in our extended high-z DLA sam¬ 
ple (INoterdaerne et al.l l2008l a nd the new ESDLA s) with val¬ 
ues in the local Galactic ISM (ISavage et al.lll97'^ . in the SMC 
(IWeltv et alJl20I2t) . and in DLAs associated with y-ray burst 
afterglows (GRB-DLAs). In the overall population, we clearly 
see a bimodality in the distribution of N(H 2 ): most detections 
have log A(H 2 ) >17, far above the typical detection limits (a few 

‘ Either because of the detec tion of Ci absorption dSrianand et all 
l200^jNoteMaeme_et^ l2010l) . the presence of 21 cm absorption 
)ISriananTeraniToi2l) . or direct evidence of H 2 lines dBalashev et al.l 

|20ll^ 


Noterdaeme_et^ 20081 : iBalashev et alJl2()TS or possi- 


Jorgenson et al.l 20141) . The increase in the fraction of 


times lO'^cm^^). In the following, we denote as “strong” (resp. 
“weak”) the systems with log A(H 2 ) > 17 (resp. < 17). The right 
panel shows the distribution of systems in each of these pop¬ 
ulations as a function of the H i column densitjQ. We hnd that 
H 2 is detected with column densities higher than 10^’ cm“^ in 
four (hve if we include the possible H 2 detection in the DLA to¬ 
wards JO 154-1-1935) ESDLAs out of seven. This is signihcantly 
higher t han the value seen in the the overall DLA po pulation 
(~ 10%, 
bly less, 

strong H 2 systems is signihcant but not as sharp as is seen in the 
Milky Way or in the Small Magellanic Cloud. In addition, the 
overall molecular fractions remain modest (~ 1% or less). 

To explain this, it must be noted that the multi-phase na¬ 
ture of the neutral gas is not equivalently probed by the differ¬ 
ent samples. The values corresponding to the Milky Way come 
from lines of sight towards nearby stars that are located only 
within ~100 pc. These should therefore probe a single cloud 
that produces most of the total observed column density and 
in which the A(Hi) and N(H 2 ) can be directly related by mi¬ 
crophysics. The situation is already different towards stars in 
the Magellanic clouds for which the observed column densi¬ 
ties may include g as from different c louds or phases along the 
same line of sight. I Weltv et aP (1201 2h also argued that previous 
A(H i) determinations in the SMC were overestimated because 
they were derived from 21 cm emission, which averages struc¬ 
tures in the ISM at scales smaller than the radio beams. Indeed, 
once the A(H i)-values are more accurately determined using Ly- 
a-absorption (i.e. along the same pencil-beam line of sight as 


^ For simplicity only the total H 2 column density is used for a given 
DLA. In spite of this, the total A(H 2 ) is dominated in all cases by a few 
components that individually also have log A(H 2 ) > 17. 
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used for H 2 measurements) higher molecular fractions are found 
in the SMC, revealing a clearer segregation between the strong 
and weak H 2 populations around logA^(Hi) ~ 21 0. In DLAs, 
the HI and H 2 column densities are measured through UV ab¬ 
sorption along the same line of sight. However, a single quasar 
sight line likely samples multiple gas components having differ¬ 
ent physical conditions, as seen from the excitation of different 
specie s (e.g. lSrianand et al.ll2005l:lLisztil2015T:lNoterdaeme et aO 
l2015h . In addition, at a given redshift, different DLAs probe dif¬ 
ferent galaxies with their own sets of physical conditions, which 
may contribute to s moothing the observat ion of any underlying 
transition. Recently. iBalashev et al.l(l2015h have used chlorine to 
show that the local metallicity and molecular fraction in the H 2 
components could be much higher than the line-of-sight aver¬ 
aged value, although this does not tell us whether the remaining 
H I is located in outer layers or in unrelated interloping clouds. 

Our results show that a large amount of Hi in ESDLAs 
could indeed be unrelated to H 2 . This is also supported by 
the similar H 2 column densities seen in several much lower 
A(Hi) systems. The large A(Hi) probably results from a low 
impact parameter of the line o f sight relative to the galactic 
centre (iNoterdaeme et alJ I20l4) where the covering factor of 
H 2 -bearing gas would be h igher owing to higher ISM pressure 
dBlitz & Rosolowskvll2006ll . 

The situation could be similar along the lines of sight to¬ 
wards afterglows of long-duration y ray burst s (GRBs) where 
DLA s are often seen with logA(Hi) > 22 (e.g. l.Takobsson et^ 
I2OO6I1. As these GR Bs are linked to the death of a massive star 
dBloom et al.lll99^ . they are probably related to star forming re¬ 
gions that are typically denser and closer to the cen tre of the 
host galaxy than quasar-DLAs dPontzen et al.l l20f(ih . Because 
GRB-DLAs may be sub ject to a very intense UV radiation field 
dTumlinson et al.ll2007h one has to exercise caution when com¬ 
paring them with quasar-DLAs. Nevertheless, although the sam¬ 
ple is still small it appears that the detection rate is consistent 
with that seen in quasar-DLAs albeit with larger molecular frac¬ 
tions. This further supports the idea that most high column den¬ 
sity lines of sight likely probe the central regions of a galaxy. 

E SDLAs are very rare and huge surveys are needed to hnd 
them dNoterdaeme et al.ll2009l 12013) . However, one could ques¬ 
tion the fact that the H i to H 2 transition may induce a bias in the 
selection of the quasars against the detection of the correspond¬ 
ing systems. 


3. The effect of ESDLAs on the colours of the 
background quasar. 

In Eig. 13 we identify the H 2 -bearing DLAs in the A(Hi)- 
metallicity plane. Interestingly three quasar-DLAs and four 
GRB-DLAs are now known beyond the limit for significant 
dust obscuration proposed by Boisse_e^^ Jl99^ and long dis¬ 
cussed in the literature (e.g. lNeeleman et al.ll2013[). Six of these 


DLAs show self-shielded H 2 . As observed byJPetitjeaii_et^ 
( 2006ll and as predicted by some models (e.g. iKrumholz et^ 
I 2 OO 9 I: ISternberg et al.l l2014t) . the H 2 detection rate is higher 
at high metallicity. Howeyer, the metallicity at which H 2 is 
found increases with decreasing A(Hi). The presence of H 2 
could be more closely related to the column density of dust 
grains (INoterdaeme et al.l l2008h : using the column density of 


^ We caution that part of the background stars were targeted for stellar 
studies (hence generally probing low extinction lines of sight), while 
others were specifically targeted to study the properties of dust, provid¬ 
ing highly reddened lines of sight. 



log N(HI) 


Fig. 2. Metallicity versus H i column density. Data samples and corre¬ 
sponding symbols are as in Fig. |3 Filled (resp. unfilled) circles repre¬ 
sent systems in which H 2 has been detected with log A(H 2 ) > 17 (resp. 
log A(H 2 ) < 17). The dotted (resp. dashed) line represents a constant 
log A(Znii) = 12.5 (resp. logAl(Znii) = 13.15). 


undepleted elements as a first-order proxy for that of dust 
(IVladilo & Perouxll2nn.5h . we can see that 10 of the 15 systems 
above a line of constant log A(Znii) = 12.5 have logA(H 2 ) > 17, 
while this fraction is only 4/66 below. 

Continuum absorption by dust and the absorption from lines 
in the Lyman series of H i and Lyman and Werner bands of H 2 
can significantly affect the quasar transmitted flux in the different 
bands when column densities become very large. We quantify 
these effects by calculating the transmission for different H i and 
H 2 column densities and different reddening. For each absorp¬ 
tion situation, the induced colour changes depend on the absorp¬ 
tion redshift, the filter responses, and the input spectrum (quasar 
continuum plus Ly a forest). For simplicity, we fixed Zabs = 2.35 
(i.e. the redshift of our strongest ESDLA towards J2140-0321), 
used the filter re sponses of the Sloan Digital Sky Survey (SDSS, 
lYork et aklbPOOl) . and considered a flat quasar spectrum. Our re¬ 
sults are shown in Fig. |3] We empirically checked that assuming 
a flat spectrum has little effect on the results. To this end, we in¬ 
troduced fake absorbers with known properties (A(Hi), A(H 2 ), 
and dust) in real non-BAL quasar spectra (with emission redshift 
close to that of J2140-0321) and derived the colour changes. We 
find very good agreement (~ 0.01 mag for r, i, z and 0.05 mag in 
g) with our simple model Q 

In the case of J2140-0321, we find that the damped Lycr 
line alone severely affects the g-band (by 0.18 mag). Similarly, 
the strong H 2 absorption lines raise the M-band magnitude by 
~0.26 mag. The importance of these line absorption is similar to 
that of the continuum absorption owing to the presence of dust 
(for E(B-V) ~ 0.05, est imated through SED profile fitting, see 
INoterdaeme et al.ll2014l) in these two bands. While the overall 
colour excesses estimated for J2140-0321 are likely not large 
enough to push the quasar out of typical colour-selections, it is 
still significant. We note that a DLA with the same characteris¬ 
tics as those of the DLA associated with GRB 080607 would in 


This consistency check is not possible in the ji-band because the cor¬ 
responding wavelength range is not covered by the SDSS spectra. 
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Fig. 3. Colours of J2 140-0321 compared to those of other non-BAL 
quasars from SDSS-111 at the same redshift. The arrows in the top right 
comer illustrate the estimates of colour excess due to Hi, H 2 , and dust 
(assuming SMC extinction law) at = 2.34. We note that because 
transmission is multiplicative, the corresponding magnitudes are ad¬ 
ditive and so are the colour vectors. The red arrow shows the esti¬ 
mated colour excess of J2 140-0321 (circled red) due to the interven¬ 
ing ESDLA. For comparison, the colour excess that would produce a 
DLA with same column densities and reddening as those associated 
with GRB080607 would reach the full ranges covered by both axes. 


turn very likely escape current colour-based selections. In addi¬ 
tion, even if such a spectrum were to be obtained, then the DLA 
would still be hard to recognise using automated detection pro¬ 
cedures because of the little transmitted flux remaining between 
damped Hi and H 2 absorption lines. This shows that alternative 
selection of quasa r s, based for example o n NIR photometry (e.g. 
iFvnbo et ^l2013l:rKrogager et al.ll2015h . are important in order 
to avoid biasing against intervening systems beyond the neutral 
to molecular transition. A promising tech nique for identifyin g 
these systems is to search for C i absorption (iLedoux et alj201^ . 


4. Conclusions 

We have extended the study of H 2 in DLAs to the very high 
Hi column density end, allowing us to uncover a significant 
increase in the fraction of strong H 2 systems (that we define 
as having logA(H 2 ) > 17) at logA(Hi) > 21.5. While the 
high A(H i)-threshold is qualitatively consistent with expecta¬ 
tions from theoretical models describing H 2 microphysics, the 
mean molecular fraction in these systems remains relatively low. 
This can be explained by the quasar lines of sight having long 
path lengths through galaxies. In this picture, most of the Hi is 
due to clouds unrelated to the molecular phase probed by H 2 . 
The threshold for local Hi to H 2 conversion in high-z DLAs 
could actually occur at A(H i) and metallicities similar to those 
in the Milky Way disc. The large H 2 column densities observed 
in EDLAs (with log A(Hi) > 21.7) could simply be due to the 
line of sight pass i ng clo ser to the galaxy centre as shown by 
iNoterdaeme et al.l (l2014l) . where the ISM pressure is higher and 
so is the probability of intercepting a molecular cloud. This is 
also consistent with the skewed A(H ^-distributions observed in 
samples of absorbers selected for their high molecular content. 


The A(H ^-distributio n of ~20 strong H 2 ab sorbers directly se¬ 
lected from the SDSS (iBalashev et al.llTOldl) is indeed biased to¬ 
wards high A(Hi) systems. Similarly. iLedoux et al.l (1201 5h ob¬ 
serve an excess of strong A(H i)-systems among C i-selected ab¬ 
sorbers, which appear to harbour high molecular content. We 
note that high m olecular content is also found in some low A(H i ) 
absorbers (e.g. ISrianand et al.l l2008t INoterdaeme et al.l 1201 Oli . 
which shows that the conversion from atomic to molecular hy- 
drogen due to micr ophysics (which occurs on pc-scales, e.g. 


Srianandet^DOlA ) does not require very high A(H i) (see also 


Muzahid et alJl2015l). 


We have investigated the impact on the quasar colours by 
the presence of systems beyond the neutral to molecular transi¬ 
tion and showed that selection of quasars with NIR photometry 
would be important in order to avoid biasing against the detec¬ 
tion of systems with high molecular content. 
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